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Abstract. Fuse spectra of the four brightest H ll regions in M 33 show absorption by interstellar gas in the Galaxy 
and in M 33. On three lines of sight molecular hydrogen in M 33 is detected. This is the hrst measurement of 
diffuse H 2 in absorption in a Local Group galaxy other than the Magellanic Glouds. A quantitative analysis is 
difficult because of the low signal to noise ratio and the systematic effects produced by having multipie objects 
in the Fuse aperture. We use the M 33 Fuse data to demonstrate in a more general manner the complexity of 
interpreting interstellar absorption line spectra towards mnlti-object background sources. We derive H 2 column 
densities of « 10*^® to 10*^^ cm~^ along 3 sight lines (NGC 588, NGC 592, NGC 595). Because of the systematic 
effects, these values most likely represent upper limits and the non-detection of H 2 towards NGG 604 does not 
exclude the existence of significant amounts of molecular gas along this sight line. 
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1. Introduction 


Observations with the Copernicus satellite have shown 
that the ISM of the Milky Way contains large amounts of 
molecular hydrogen. The mass in H 2 may even dominate 
that in atomic form (see Savage et al. 1977). Since 1996 
other galaxies also have become accessible for investiga¬ 
tion of the content of H 2 , using Orfeus and Fuse far ul¬ 
traviolet spectra. Both the Large and the Small Magellanic 
Cloud appear to have substantial amounts o f H 2 as found 
from absor ption lin e stud ies (de Boer et al. 1998 ; Richter 
et al. 1998 Richter ^000 , Tumlinson et al. |2002 


Fuse has sufficient sensitivity to record spectra of 
bright sources in other nearby galaxies such as M 31 and 
M 33. It is of importance to extend our knowledge about 
the content in molecular gas to galaxies other than the 
Milky Way and the LMC and SMC. These galaxies have 
different histories, different stellar contents and different 
abundances. An intercomparison of these galaxies may 
lead to clues about the interplay of molecular cloud and 
star formation with chemical content and radiation field. 


M 33 is the smallest spiral galaxy of the Local Group at 
a distance of ~ 800 kpc. It has a ragged optical appearance 
and bright Hii regions, some of them catalogued as NGC 
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* Based on observations made with the NASA-CNES-CSA 
Far Ultraviolet Spectroscopic Explorer, available in the public 
archive. FUSE is operated for NASA by the Johns Hopkins 
University under NASA contract NAS5-32985. 


objects. The galaxy is inclined at 55° and shows signs of 
bending as derived f rom 2 1 cm radio synthesis observa¬ 
tions (Rogstad et al. 1976 ). In the direction of two of the 
bright H ii regio ns CO has been found in emission (Wilson 
& Scoville 1992 ). The metallicity of M33 is subsolar with 
a radial abundance gradient which differs for different ele¬ 
ments. For example, neon abundances in the observed H ii 
regions vary bet ween 0.2 and 0.3 dex subsolar (see Willner 
& Nelson-Patel 2002|) . Vibrationally excited H 2 has been 
found in emission by Israel et al. (1990) towards NGC 604, 
most likely from background gas. 

In this paper we investigate Fuse spectra of four bright 
objects in M 33. Each of these consist of a number of young 
and thus UV-bright stars surrounded by an H ii region. In 
3 of the 4 spectra H 2 is detected in absorption in M 33. 

In Sect. 2 we give some information about the instru¬ 
ment, the targets, and the data reduction, in Sect. 3 the 
difficulties with the interpretation of interstellar absorp¬ 
tion spectra towards multi-object background sources are 
described. The measurements are presented in Sect. 4, fol¬ 
lowed by a discussion of the results in Sect. 5. 


2. Data 

2.1. The instrument 

Fuse is equipped with four coaligned telescopes and 
Rowland spectrographs. The detectors are two microchan- 
nel plates. For a description of the instrument and its on- 
orbit performance see Moos et al. (2000D and Sahnow et al. 
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Fig. 1. The Fuse 30 x 30 arcsec LWRS aperture overlaid on the four Hii regions / star clusters. For NGC 588, 
NGC 595, and NGC 604 HST WFPC 2 images (central wavelength 1730 A) are shown, for NGC 592 only a DSSl 
image was available. 


(2000). Information about Fuse data and their analysis 
can be found in the Fuse Data Handbook]^. 

The typical velocity resolution of Fuse spectra is 
about 20 kms“^. Because of variations in the alignment of 
the 8 different detector channels the measurement of abso¬ 
lute radial velocities is problematic and requires a separate 
discussion (see Sect. 4.2). 

The sensitivity threshold of Fuse is set by the un¬ 
certainty in the background intensity level. It is below 
(2 — 4) X 10“^^ ergcm“^ s“^ A“^. To allow a reasonable 


Available at ittp://fuse.pha.jhu.edu/analysis/dhbook.htm' 


short integration time the source flux should be higher 
than about 1 x 10“^^ ergcm“^ s“^ A“^. 


2.2. The targets 

The targets are the UV-bright cores of the Hii regions 
NGC 588, NGC 592, NGC 595, and NGC 604 in M 33. For 
the Fuse observations an entrance aperture with size of 
30 x 30 square arcsec was used, being equivalent to 120 x 
120 pc^ in M 33. The precise locations of the apertures 
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Fig. 2. For each line of sight a sample of five H 2 and four metal lines is shown. Because of peculiarities in the spectra 
(sight line structure, noise features, structures in the background source spectrum) the choice of transitions is not the 
same for all four sight lines. The dotted lines mark the radial velocities of Galactic and M 33 gas. Identified other 
transitions within the range of the plots are labeled in the spectra as follows: H 2 Lyman: Oa: R(0),4-0; Ob: R(0),6-0; 
la: R(l),l-0; lb: R(l),2-0; Ic P(l), 3-0; Id: R(l),4-0; le: R(l),5-0; 2a: R(2),2-0; 2b: P(2),3-0; 2c: R(2),4-0; 2d: 
R(2),5-0; 2e: P(2),5-0; 3a: R(3),4-0; 3b: P(3),6-0; 4a: R(4),3-0; 4b: R(4),4-0; 4c; P(4),5-0; 5a: P(5),4-0; 5b: R(5),5-0; 
5c: P(5)(5-0). Metals: ma: An (1048 A); mb: Fell (1055 A); me: Feiii (1122 A); md: Sill (1020 A). M 33 components 
are marked with a star. 
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Table 1. Data on Fuse pointings on M 33 Hii regions. The position angle is the angle of the Fuse y axis (east of 
north). N 20 is the number of stars having an UV flux (at 1730 A) of at least 20% of the brightest source in the area 
covered by the Fuse aperture. For NGC 592 A ^20 is unknown because no public HST image is available (see Fig. |^). 


Object 

Observation ID 

02000.0 

<52000.0 

Pos. angle 

Exp. time 

N 20 

NGC 588 

a0860404001 

a0860404002 

0132 45.50 

-b30 38 55.00 

334.9° 

2819 s 

2401 s 

8 

NGC 592 

a0860202001 

a0860202002 

0133 12.30 

-b30 38 49.00 

334.9° 

2164 s 

1800 s 

7 

NGC 595 

a0860303001 

0133 33.60 

-b30 41 32.00 

322.4° 

3620 s 

12 

NGC 604 

a0860101001 

a0860101002 

0134 32.50 

-b30 47 04.00 

323.1° 

3633 s 

3520 s 

9 



log (Na) log (Na) 


Fig. 3. Curves of growth for absorption by H 2 in M 33. The solid curves are for the best fit (or best guess) 6-value, 
the position of the data points on the log(/A) axis corresponds to the “optimum” column densities given in Table |[ 
The dashed and dotted curves represent the upper and lower limits of 6-values consistent with the data. They were 
used to determine the minimum and maximum column densities given in Table |[ 


with their orientation and the exposure times are given in 
Tabled 

Inspection of imagery of the M 33 H ii regions shows 
that the Fuse aperture has included numerous objects, of 
which only few were bright enough in the UV to have con¬ 
tributed to the background flux. Nevertheless, the signal 
recorded is the sum of the flux of several stars. This implies 
that the spectra provide a blend of spectra from different 
lines of sight with different line of sight gas (composition 
and velocity) structure, in particular within M 33. 

In an attempt to account for this complexity, the im¬ 
agery (Fig. 0 ) has been investigated to identify the bright¬ 


est sources in the field. We have included the scanty in¬ 
formation from that in Table |^. 

The complexity of the background sources complicates 
the interpretation of the interstellar absorption features as 
described in Sect. 3. Furthermore, the shape of the con- 
tinua is difficult to assess in some cases. 

While the FUV fluxes of all 4 H ii regions 
are well above the sensitivity of Fuse (« 5 x 
10“^'* ergcm“^ s“^ for NGC 588, NGC 592, and 
NGC 595), the exposure times were too short to have 
more than 20 to 50 counts per 20 kms“^ resolution el¬ 
ement. An exception is the NGC 604 spectrum where the 
flux (« f X 10“^^ ergcm“^ s“^ A“^) provides roughly 400 
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Fig. 4. Curves of growth for absorption by metals in M 33. Curves for the adopted 6-value (solid lines), as well as for 
the upper and lower limits (dashed and dotted lines), are shown. 


counts per resolution element. Sample interstellar spec¬ 
tra are shown in Fig. |^, a few other spectra are given by 
Wakker et al. ( ^002] ). Note that the spectra in Fig. ^ are re¬ 
binned over 10 pixels. The resulting pixel size corresponds 
roughly to one resolution element. 

2.3. Data reduction 

The raw spectra were calibrated with the Calfuse 1.8.7 re¬ 
duction pipeline. For analysis the spectra were smoothed 
with a 5 pixel (NGC 604) or a 10 pixel boxcar fil¬ 
ter (NGC 588, 592, 595) in order to handle the strong 
noise. Equivalent widths were measured using trapezium 
or Gaussian fits. The errors were calculated taking into 
account the photon statistics and the uncertainty in the 
choice of the continuum. 


A calibration of the velocity zero point was made by 
comparison to radio data (see Sect. 4.2). 

3. General considerations for the interpretation 

3.1. Velocity structure 

Due to the presence of several background targets within 
the Fuse aperture the resulting spectrum is the sum of 
single spectra with different fluxes and absorption fea¬ 
tures. At the distance of M 33 one arcsecond corresponds 
to a projected distance of about 4 pc, so the Fuse aper¬ 
ture covers an area of roughly 120 x 120 pc^. As Fig. 
shows, the stars which dominate the UV fluxes lie within 
smaller areas of about 50 x 50 to 70 x 70 pc^. On this scale 
significant fluctuations in velocity structure, 6-values, and 
column densities of the M 33 interstellar medium can be 
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Fig. 5. The Fuse spectra contain blends of the spectra 
of the brightest objects contained in the aperture. The 
effects of such blends are explored assuming two equally 
bright background stars but one having a 2.6 dex higher 
H 2 column density. Panel A) shows the curve of growth 
for the single sight-line spectrum with low H 2 column den¬ 
sity, panel B) the one for high H 2 column density. Panel 
C) shows the optical depth relation for the added spec¬ 
trum. In each panel the column densities of the rotational 
levels are indicated in brackets, in C) the mean values 
of the single sight line column densities are given. It is 
obvious that in panel C) a single cloud curve of growth 
interpretation is inappropriate. 


expected, especially for molecular gas that exists in con¬ 
fined regions. 

The velocity structure seen in such a sum spectrum 
thus reflects the combination of the possible presence of 
gas clouds with different velocities in front of the clus¬ 
ter, and/or possibly differing absorption velocities on the 
individual sight lines. 

For the Galactic foreground absorption the latter effect 
of widely separated background sources in M 33 is proba¬ 
bly of minor importance. While column density variations 
in the ISM are known to exist on small scales, velocity 


variations appear to be rather small (see e.g. Andrews et 
al. ^OOll) . 

Background sources separated in angle have yet an¬ 
other effect on the velocity structure of the spectrum. 
If the background sources are spread in the direction 
of dispersion, equal intrinsic absorption velocities and 
structures would be spread out to a maximum of about 
100 kms“^. As can be seen in Fig.|^ the H 2 and metal line 
absorption is wider than normally seen, but substantially 
less than the indicated maximum spread. We conclude 
that the velocity smearing due to multiple background 
sources is less than Ri 50 kms“^. 


3.2. Absorption strength 

In order to judge in which way the blending of the light 
from the individual objects affects the interpretation of 
the spectra we will consider a simple but realistic case. 

For a single line of sight, several clouds and velocity 
components contribute to the interstellar absorption spec¬ 
trum and the optical depth in an absorption profile is the 
sum of optical depths along the sight line. 

In contrast, for a sum of spectra of different sight lines 
the measured equivalent widths of interstellar lines are 
mean values of the equivalent widths belonging to indi¬ 
vidual sight lines weighted by the flux of the background 
sources. 

We have constructed a simple case of blending to in¬ 
vestigate in which way the data become less accessible 
for interpretation. We consider two sight lines which con¬ 
tribute to one total spectrum by equal amounts. The ab¬ 
sorption structure of the H 2 lines on these lines of sight 
has been chosen to be different, one has a H 2 column den¬ 
sity of 4.4 X 10^® cm“^, the other has 1.8 x 10^® cm“^. 
Consequently, the lines in one spectrum are much weaker 
(if not even ‘absent’), and the light of that star blends 
some of the absorption present in the other spectrum. 
Clearly, the equivalent widths derived from the sum spec¬ 
trum have been altered compared to those of the single 
spectra. 

Column densities determined from lines in the linear 
part of the curve of growth are the weighted mean of the 
column densities on the individual sight-lines. If lines in 
the damping part of the curve of growth are added, the 
measured column density is proportional to 

( 1 ) 


This is less or equal to 
1 


Si 9 i 




2 

i f 


( 2 ) 


the weighted mean of the individual column densities, 
where Wi are the equivalent widths in the individual spec¬ 
tra and gi are the weights. In the logarithmic part of the 
curve, where the equivalent widths depend strongly on the 
&-value the measured column densities are again another 
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Table 2. Column densities logiV {N in cm and velocities of M 33 absorption components detected in the Fuse 
spectra. Note that the H 2 column densities are “best guess” values based on the “optimum” 6-values given in Table |[ 



NGC 588 

NGC 592 


NGC 595 

NGC 604 




A 

B 


A 

B 

r’LSR [kms 

-140 

-150 

-110 

-160 

-235 

-140 

Fell 


14.951°:?® 

15.001°:®° 

15.101°:?° 

i5.ool?:“ 

14.601?:“ 

An 



1 9^+2.45 

10.ZO_Q 5Q 

14.401°:?® 

13.65l?:?®0 

< 12.9 

Oi 





16.20l?:“ 


H 2 

16.6 

17.4 

- 

17.3 

< 15.0 


kind of mean value. The situation is even more compli¬ 
cated if line profiles from different parts of the curve of 
growth are mixed. 

We illustrate the effects of the considered case in Fig. ||. 
From this simple example, considering the resulting opti¬ 
cal depth relation (see Fig. || c), we conclude that column 
densities can be derived from the Fuse M 33 H ii region 
spectra only in an approximate manner. In particular, a 
single 6-value cannot be used to interpret the absorption 
line strengths (as visible in Fig. ^ c). 


3.3. Extinction 


If the dust abundance in front of a cluster is spatially inho- 
mogenous, a selection effect has to be considered for the 
interpretation of measured column densities. In general 
the dust abundance is correlated with the column density 
of the gaseous ISM. Sight lines with higher gas column 
densities have a higher colour-excess, reducing the far UV 
flux in the spectra and thus lowering the weight in the 
sum of spectra. Therefore, column densities measured in 
the composite spectrum of a cluster can be expected to 
be lower than the actual (un-weighted) average column 
densities towards the cluster. Because the abun dance of 
H 2 is correlated with dust ( see Savage et al. 1977, Richter 


200 c , Tumlinson et al. 2002| ), sight lines with high H 2 col¬ 
umn densities are also expected to be suppressed in the 
composite spectrum. 


4. Results 

4.1. Column densities 

4.1.1. Molecular hydrogen 

Despite the systematic errors described in Sect. 3 we use 
single cloud curves of growth for column density estimates 
because the present data do not allow a more detailed 
modelling. For each rotational level the measured equiv¬ 
alent widths form a set of data points with unknown po¬ 
sition on the log{NfX) axis. These sets are shifted along 
this axis in order to fit a single-cloud curve of growth. 

In the present case the range in log(/A) covered by 
each set of data points is small and the uncertainties in 
the equivalent widths are large. Thus the ranges of possible 
6-values and column densities are substantial. 


As mentioned in Sect. 3.2, in the case of a multi¬ 
ple background source, column densities derived from the 
damping part are systematically lower than the weighted 
mean of column densities for the individual sight lines. 
Thus the upper limits are not as rigid as they would be 
for a single sight line. 


4.1.2. Metals 


The curve of growth method was also used to obtain Fe ii 
and An column densities. The results, especially for Fell, 
can be assumed to be more reliable than those for H 2 . 
Due to the usually higher 6 - values in the warm neutral 
medium and the large range of log/A - values of acces¬ 
sible Fell lines some data points lie on the linear part of 
the curve of growth (see Fig. ^). Fell and An are the 
predominant ionization states of Fe and Ar in the neu¬ 
tral interstellar medium. The cosmic abundance ratio is 


log(Fe/Ar)~ 1.1. According to Sofia & Jenkins (199S) ar¬ 
gon is unlikely to be bound to dust in clouds that are thin 
enough to be observed in absorption. In partially ionized 
clouds, significant amounts of argon can be ionized due 
to its large ionization cross section although its ionization 


energy lies above that of hydrogen. Sofia & Jenkins (1998) 


measured An depletions between 0.2 and 0.6 on Galactic 
sight lines. Therefore, in general, the Fe/Ar ratio will yield 
an estimate for the iron depletion that is too low. 

The degree of ionization of oxygen is coupled to that 
of hydrogen due to charge exchange reactions. Thus it 
would be desirable to know O i column densities, which 
then could be converted into H i column densities provided 
the metallicity is known. Unfortunately the 1039 A line 
of OI is affected by blending with H 2 lines and emission 
features and lies in the flat part of the curve of growth. 
We have included a plot of this line in Fig. ||, but do not 
analyse it. Only the spectrum of NGC 604 has a sufficient 
count rate below 1000 A for a measurement of the O i line 
at 976 A. 

Metal column densities are given in Table |^. 


4.2. Velocities 

Because of the possible apparent velocity shifts induced 
by multiple background sources, and the variations in the 
alignment of the optical channels of Fuse we measured 
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only velocity differences between Galactic and M 33 ab¬ 
sorption components. Noise and continuum uncertainties 
limit the precision to about 10 kms“^. 

Using data from the Leiden - Dwingeloo - Survey 
(Hartmann & Burton 1997 ), the velocity of the Galactic 
Hi in the direction of M 33 can be determined. We fit 
fLSR = 0 kms“^ and convert the measured Fuse veloci¬ 
ties to LSR velocities. 

We compared emission line data (CO, Hi, Hii or other 
ions) from the literature with the Fuse absorption line 
data to get hints on the location of the absorbing gas rel¬ 
ative to the H II regions and the M 33 disk. Velo cities of H i 
21 cm emission are taken from Rogstad et al. ( 1976 ) who 
found two disk components, a “main” and a more patchy 
“weak” one. These Hi “main component” velocities are 
in agreement with velocities extracted from the work of 
Deul & van der Hulst ( |1987 ) based on observations with a 
5 to 10 times higher angular resolution. Notably, no “weak 
component” shows up in the Deul & van der Hulst data 
at the positions of the H ii regions. This may be due the 
patchiness of that component or to lower sensitivity of the 
higher resolution data. 


4.3. Comments on the sight lines 
4.3.1. NGC 588 

One M 33 velocity component is detected in absorption 
by metals and H 2 in this direction. However further com¬ 
ponents could be hidden due to the limited resolution and 
the noise. The radial velocity of —140 kms“^ is about the 
same as that of the weak H i component seen in emission. 
The velocity of the main Hi component (—160 kms“^) is 
more similar to that of the emitting g as in t he H ii region 
(—168 kms“^, see O’Dell & Townsley 1988 ). The column 
densities of Fe ii and Ar i indicate some depletion of iron 
and thus the presence of dust. 


4.3.2. NGC 592 

For this Hii region no emission line results are avail¬ 
able. The smoothed Hi velocity maps give about —160 
and —146 kms“^ for the “main” and the “weak” com¬ 
ponent, respectively. The Fuse absorption velocities of 
— 150 and —110 kms“^ are closer to the velocity of the 
“weak component”, but —150 kms“^ is also consistent 
with the main H i component velocity. H 2 is found only in 
the —150 kms“^ component. 


4.3.3. NGC 595 

An absorption component at w —160 kms“^ is detected 
in metal and H 2 lines. Absorption lines by higher rota¬ 
tional levels (J = 2, 3,4) are possibly shifted to more neg¬ 
ative velocities by 10 or 20 kms“^. Despite the rather 
high “best guess” value for the H 2 column density, the 
Fell to An ratio points to mild depletion and relatively 
low dust abundance. CO emission observations by Wilson 


Table 3. M 33 H 2 column densities [cm“^]. For none 
of the sight lines a definite 6-value could be determined. 
Thus column densities of the rotational levels for maxi¬ 
mum,optimum,and minimum 6-values [kms“^] consistent 
with the data are given in the 3rd, 4th, and 5th column. 
For NGC 588 the damping limit was used to derive the 
upper limits of the column densities. Though the actual 6- 
value is certainly > 0 kms“^, a vanishing Doppler param¬ 
eter cannot be excluded by the data. For NGC 604 upper 
limits were derived using the Doppler- and the damping- 
part of the curve of growth. 



J 


log A 




^max 

^opt 

^min 

NGC 588 


23.0 kms“^ 

10.0 kms"^ 

0.0 kms“'^ 


0 

15.6 

16.1 

18.3 


1 

15.7 

16.4 

18.5 


2 

15.2 

15.6 

18.1 


3 

14.7 

15.0 

17.9 


4 

14.6 

14.9 

17.8 


5 

< 14.5 

< 14.6 

< 16.8 


6 

< 14.5 

< 14.5 

< 16.7 


7 

< 13.7 

< 13.7 

< 15.2 

NGC 592 


17.0 kms"^ 

10.0 kms"^ 

4.0 kms“^ 


0 

16.0 

17.1 

18.4 


1 

15.8 

16.9 

18.4 


2 

15.6 

16.4 

18.2 


3 

15.1 

15.4 

17.2 


4 

15.0 

15.2 

17.1 


5 

< 14.9 

< 15.0 

< 15.5 


6 

< 14.6 

< 14.7 

< 15.2 


7 

< 14.7 

< 14.8 

< 15.3 

NGC 595 


12.0 kms”'^ 

4.5 kms“^ 

1.5 kms“^ 


0 

15.5 

17.0 

17.9 


1 

15.7 

16.9 

18.1 


2 

15.2 

16.5 

17.7 


3 

14.8 

15.1 

17.1 


4 

14.7 

14.9 

16.6 


5 

< 14.8 

< 15.1 

< 17.1 


6 

< 14.9 

< 15.5 

< 17.3 


7 

< 14.6 

< 14.8 

< 16.7 

NGC 604 


inf 


0 kms“^ 


0 

< 13.7 


< 15.7 


1 

< 14.3 


< 15.9 


2 

< 14.4 


< 16.7 


3 

< 14.0 


< 15.8 


4 

< 14.4 


< 16.6 


5 

< 14.0 


< 15.7 


6 

< 13.8 


< 15.3 


7 

< 13.8 


< 15.5 


& Scoville ( 1992| ) revealed the existence of two molecu¬ 
lar clouds (vlsr = —186 kms“^ and —189 kms“^) south 
of the cluster but within the region covered by the Fuse 
aperture. However, clouds detectable in CO emission nor- 
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mally have so much dust that UV star light from behind 
should not get through. So our data do not refer to such 
dense and opaque gas. 

The velocity of the main H i emission component is 
~ —183 kms“^, that of the “weak” Hi component is 
w —155 kms“^. Possibly the w —155 kms“^ component 
dominates in H 2 and metals while in the w —183 kms“^ 
component the higher J-levels of H 2 contribute more (see, 
e.g., the similar case of HD 5980 in the SMC; Richter et 
al. [I9^. 


4.3.4. NGC 604 

Two absorption components are found in metal lines: 
a strong one at w —235 kms“^ and a weaker one at 
w —140 kms“^. The velocity of the first one is consis¬ 
tent with the velocity of the main H i emission component 
(u « —236 kms“^), while the velocities of the weaker ab¬ 
sorption component and the weak H i emission component 
(u « —206 kms“^) do not match. 

No absorption by H 2 at M 33 velocities is present, al¬ 


though Wilson & Scoville (1992) detected CO emission 
south-east of the area covered by the Fuse aperture at 
a velocity of —242 kms“^. The LSR velocity of the Hii 
region is —246 kms“^ (O’Dell & Townsley 1988), similar 
to that of CO shifted about -1-10 kms“^ compared to the 
main Hi emission component. 

The comparatively high count rate in this spectrum 
allows the determination of a relatively rigid upper limit 
for the H 2 column density. The non-detection of H 2 in 
absorption is consistent with the measured Fell and An 
column densities (see Table |^. In both velocity compo¬ 
nents the abundances do not indicate any depletion. For 
component A it was possible to measure the O i line at 
976 A. The O i to Fe ii ratio is consistent with a cosmic 
abundance ratio and no depletion of iron. 


5. Discussion 

5.1. H 2 or no H 2 

The results for NGC 588, NGC 592, and NGC 595 
are very similar: almost identical Fe ii column densities, 
log(Feii/(Ari) between 0.3 and 0.7, indicating moderate 
iron depletion, and radial velocities consistent with the 
“weak” H i component and differing from the main disk 
component. 

The velocity difference between the H 2 detected in ab¬ 
sorption and the emission from the H ii regions and the 
molecular clouds in their vicinity indicates that the ob¬ 
served diffuse H 2 is probably not located near the Hii 
regions. Towards NGC 604 the velocity of the main ab¬ 
sorption component is consistent with the velocity of the 
main H i emission component. While the iron column den¬ 
sity is about the same as towards the other targets, no H 2 
is found on this sight line. 

For NGC 595 and NGC 604 photometric results have 


been published. Wilson & Matthews (1995) give E{B — 


V) = 0.35 for NGC 604 and E{B - V) = 0.25 for 
NGC 595, based on Palomar 60 inch data. They state 
that the colour magnitude diagrams suggest the presence 
of “significantly more differential reddening in NGC 604 
than in NGC 595”. From HST observations of NGC 595 
Malumuth et al. ( |l9^ find E{B -V)= 0.36 ± 0.28. We 
conclude that the average dust abundance towards both 
clusters is comparable, with a considerable scatter which 
is probably larger towards NGC 604. Thus the extinction- 
effect discussed in Sect. 3.3 is heavily influencing the ap¬ 
parent H 2 column densities in front of the clusters. 

As the G alacti c and Magellani c H 2 surveys show 
(Savage et al. 1977 , Tumlinson et al. 2002), on sight lines 
with an E{B — V) of 0.2 or larger H 2 column densities 
between 10^® and 10^^ cm“^ are common. Though most 
of the E{B — V) = 0.36 towards NGC 595 is supposed 
to originate in M 33 ( Galactic E{B — V) k. 0.06, see 
Burstein & Heiles 1982 ), the estimated H 2 column den¬ 
sity (« 2 X 10^^ cm“^) is rather low. There is no sign of 
lines from J = 0 or J = 1 lying in the damping part of 
the curve of growth as it would be expected for high col¬ 
umn densities. A rough estimate of the E{B — V) of the 
FUV brightest stars in NGC 595 b ased on the results pub¬ 
lished by Malumuth et al. ( 1996 ) yields a value of about 
0.1. For a colour excess like this an H 2 column density of 
10^’’ cm“^ is realistic. 


The photometric data for NGC 604 are not that de¬ 
tailed as for NGC 595 but it is likely that the extinction 
selection effect is even more significant. If the scatter in 
extinction is indeed larger towards NGC 604, it is prob¬ 
able that there are some almost unreddened stars which 
outshine the rest of the cluster in the far UV. Independent 
of the actual average column density of H 2 in front of the 
cluster, the few sight lines with negligible amounts of dust 
or H 2 will dominate the UV spectrum and thus lead to a 
non-detection of H 2 . 

Our results make clear that it is difficult to get precise 
quantitative information from absorption spectroscopy of 
extended or multiple background sources using a large 
spectrograph aperture. 


5.2. The bottom line 

For the first time diffuse molecular hydrogen has been 
detected in absorption in a local group galaxy farther away 
than the Magellanic Clouds. Since star formation is going 
on in M 33 and molecular clouds had been detected in 
emission by CO and excited H 2 , the presence of diffuse 
H 2 is not unexpected. The estimated column densities are 
similar to those found for Milky Way and LMC along low- 
extinction sight lines. 

More reliable data of other sight lines are desirable 
to investigate the physical conditions in the diffuse inter¬ 
stellar medium of M 33, especially parameters like radi¬ 
ation field and H 2 formation rate. So far, 10 bright OB 
stars in M 33 have been observed with Fuse. Maybe in 
at least some of these lines of sight fewer stars lie in the 
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field of view than in the data analysed here. Also, fu¬ 
ture work with the smaller MDRS or HIRS apertures of 
Fuse may alleviate the extinction-selection problem en¬ 
countered with the current data. 
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